arXiv:hep-ph/9706530v2 7 Jul 1997 


SLAC-PUB-7565 
hep-ph/9706530 
June 1997 


PION PRODUCTION FROM BAKED-ALASKA 
DISORIENTED CHIRAL CONDENSATE [] 

G. Amelino-Camelia, J. D. Bjorken [], and S. E. Larsson 

Theoretical Physics, University of Oxford 
1 Keble Road, Oxford 0X1 3NP 
United Kingdom 


Abstract 

We study the various stages of the evolution of chiral condensates disoriented via 
the “baked-alaska” mechanism, in which the condensates are described as the products 
of external sources localized on the light cone. Our analysis is based on the classical 
equations of motion of either the linear or the nonlinear sigma model. We use the asso¬ 
ciated framework of coherent states and, especially, their source functions to make the 
connection to the distribution functions for the produced particles. We also compare 
our classical approach with a mean-field calculation which includes a certain class of 
quantum corrections. 
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1 Introduction 


Recently, in order to interpret events with a deficit or excess of nentral pions observed in 
cosmic ray experiments, there has been increased interest in the conjectnre that it might 
be possible to prodnce in high-energy collisions disoriented chiral condensates (DCCs), i.e. 
correlated regions of space-time wherein the qnark condensate, {0\qL(lR\0), is chirally rotated 
from its nsnal orientation in isospin space. 

On the theoretical side there has been great interest (see, e.g., Refs. [||]-[^) both in the 
development of technical tools snitable for the description of this possible phenomenon, and 
in the exploration of the possibilities opened by DCCs as probes of the strnctnre of qnantnm 
chromodynamics, most notably in relation to the chiral phase transition. 

The idea that snch DCCs might be produced in high energy collisions at existing or 
planned hadron or heavy-ion accelerators has generated some experimental interest. In 
particular, one of us is co-spokesman for a Fermilab experiment looking for DCCs in 
hadron-hadron collisions. In high energy p-p collisions which lead to a sizable multiplicity of 
produced particles, but not necessarily with high-p^ jets in the final state, the time evolution 
is quasi-macroscopic, because the hadronization time can be rather large, up to 3-5 fm. At 
times t before hadronization, the initial state partons, produced in a volume much smaller 
than a cubic fermi, may stream outward at essentially the speed of light in all directions, 
occupying the surface of a sphere of radius t (in units such that the speed of light is 1). Most 
of the outgoing energy/momentum is expected to be concentrated near the light cone, i.e. 
on the surface of the expanding ‘fireball’. 

However, the interior of the fireball is also an interesting place. If its energy density is 
low enough, the interior should look very similar to the vacuum, with an associated non¬ 
vanishing quark condensate. Since the energy density from the intrinsic chiral symmetry 
breaking is small |^, |l^, and since the fireball surface isolates the interior from the exterior 
of the light cone, it is reasonable to consider the possibility that well inside the light cone 
the quark condensate might be chirally rotated from its usual orientation. At late times 
this disoriented vacuum would relax back to ordinary vacuum, via radiation of its collective 
modes, the pions. The properties of the radiated pions should be strongly affected by the 
semiclassical, coherent nature of the process. In particular, one may expect anomalously 
large event-by-event fluctuations in the ratio of the number of charged pions to neutral pions 
produced. Assuming that the event-by-event deviation of the quark condensate from its 
usual orientation be random, one finds []^ ||, |^, ||, |^, that the distribution P{f) of the 
neutral fraction 
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is given by 
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at large Ntot- Most notably, this implies that for 
extreme values of / is very different from ordinary pion production (which is given by a 
binomial distribution), in which the fluctuations are expected to be peaked at / = l/3 and 
fall exponentially away from the pea k. 

Experimental DCC searches |l^ are thus far largely based on the structure of Eq. (|^). 
However, it is probable that this robust property will not be sufficient on its own for the 
identification of phenomena involving DCCs. In particular, it appears that an understanding 
of the geometry of the phenomena would be very useful for experimental searches. 

In this paper we develop a description of “baked-alaska” (hot fireball surface with cooler 
and disoriented inside) DCCs which is well suited for the study of realistic geometries. Our 
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analysis is based on the classical eqnations of motion of either the linear or the nonlinear 
sigma model and we nse the associated framework of coherent states to make the connection 
to the distribntion fnnctions for the prodnced particles. 

In the next section we review the coherent-state formalism with emphasis on the associ¬ 
ated particle flnx, which is seen as directly related to classical sonrces localized on the light 
cone. Sec. 3 is devoted to the sigma model description of “baked-alaska DCCs”, dehned 
via the above-mentioned mechanism involving a hot hreball shell and its cooler (and some¬ 
times disoriented) interior. In Sec. 4 we derive a simple solntion of the classical nonlinear 
sigma model, and follow the evolntion of the associated coherent state from beginning to 
end, inclnding a derivation of the flnx of pions. In Sec. 5 we discnss how onr analysis can 
be generalized to the linear sigma model, and establish connections with a somewhat similar 
model stndied by Lampert, Dawson, and Cooper |T^. Sec. 6 is devoted to closing remarks. 


2 Particle production from classical sources 

2.1 Coherent state formalism 

In the following sections we will be investigating a problem which in the classical limit is 
described by the held eqnation 


(D-\-m^)((){x) = J{x) , (3) 

for appropriate bonndary conditions. The properties of the qnanta associated with a held, 
snch as 0 in (H), that are prodnced by a classical cnrrent sonrce J can be stndied nsing the 
coherent state dehned by 


I J) = exp 
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where J{k) is the Fonrier transform of the sonrce J{x), is the nsnal creation operator, 
and the integral measnre is given by 


dk = 


d^k 


2'7i6{k'^ — m^)0{k^) 


( 5 ) 


Note that, becanse of the integration measnre, the actnal contribntion of J{k) to the integrals 
in (^) comes from the mass-shell with kh' > 0, and sometimes it proves convenient to introdnce 
the notation 


J{k) = , (6) 

where 

ujk = \Jk‘^ + 171^ > 0 . (7) 

We remind the reader that the relation between the coherent state \J) and the sonrce 
J can be derived from the familiar (text-book) analysis of the S-matrix associated with J. 
Provided that J is snfhciently well localized in space and time (so that the idea of a scattering 
process is well dehned), this S'-matrix is given by 

S = : exp J d^y (i)f{y)J{y) : , ( 8 ) 
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where 4>f{y) is the free scalar held operator. Taking into account the normal ordering, one 
can rewrite S as 


S = exp 

If Xk \J{k)\^ 

exp 

= exp 

'iJdkimf 

exp 


dS J{y)^f{y) exp 
dk J{k)a){k) exp 


I A J{y)^~f{y) 

dk J{k)d{k) 


( 9 ) 


In particular, this implies that the “out” state 5“^|0) corresponding to a scattering process 
having the vacuum |0) as “in” state, is given by the coherent state of Eq. (^). 

The most appropriate tool for the description of the particle production associated with 
the “out” state is the generating functional, which in the present case takes the form 


G[z4,] 


exp 


dk\J{k)\\z^{k)-l) , 


( 10 ) 


from which one can obtain all the inclusive and exclusive factorial moments. In particular, 
the inclusive spectrum of particles as a function of momentum is given by 

, = . ( 11 ) 

Note that all higher correlation functions vanish for a coherent state, as is easily shown 
by taking further derivatives of the generating functional. Thus, we see that the main object 

of interest in this description is the function J{k). 


O dN^_ 
d'^k 


dG[z^] 
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2.2 Extracting J from the field 

Since one’s intuition is that the flux of particles depends directly on the held, it is instructive 
to see how the value of J{k) can be derived from the (classical) held at late times. Assuming 
that the “in” state has vanishing classical held (it is the classical “vacuum” state), (j)in{x) = 0, 
we can write 

4>out{x) = J d^ dG)[x- y) J{y) , (12) 

where D^) = Dret — Dadv is the diherence between the retarded and advanced Green func¬ 
tions. Taking the Fourier transform (and using the convolution theorem) we get 

~^out{k) = bG)[k)J{k) 

= 271 S{k‘^ — m'^) e{k^) i J{k) , (13) 

where we have substituted for the explicit form of DG^{k) in the last line. It is then useful to 
consider the three dimensional Fourier transform of (j)out{x)- We use the following conventions 

f{k) = jd^xb^^fix), 

P\t,k) = j d^xe-G^-^f{kx) , 
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which imply 
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From (^) and (pIBD it follows that the three dimensional Fonrier transform of (poutix) is given 
by 


k) = 


dk^ .-ikH 


e 27iS{k^ — m^)e(fc°) iJ{k) 
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where we have nsed the fact that J{x) is real 

J{ky = J{-k) . 

Upon observing that ([T7|) also implies 

= -i[[iJ{k)]e-^^*=^ - [iJ{k)yy^*=^] , 

we can solve for J{k) in terms of and its time derivative 


(17) 
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This relation is extremely nsefnl for problems set up in such a way that what is known is 
encoded in the equations of motion plus the initial held conhguration; in fact, (|2(]| ) allows 
us to derive an associated current source (which in turn describes the particle hux via (|^) 
from the late-time held. 


3 Sigma-model description of baked-alaska DCCs 

The 0(4) sigma-model is typically used as a model of hadron dynamics in DCC studies. It 
is simple enough to be treatable, has the correct chiral symmetry properties, and describes 
the low energy phenomenology of pions. However it must be kept in mind that it is, at best, 
a crude approximation to the chiral ehective low-energy Lagrangian of QCD. 

The Lagrangian density of the linear sigma-model is 

^ - flf + Ha . (21) 

Ff = 0 in the chiral limit, = 

A meaningful sigma-model description can start at some small proper time, of order 0.2- 
0.3 fm, near the light cone, when the collective coordinates a and vr become relevant |T^, |T^. 
At this early proper time the distribution of the chiral held 


^ = a + in ■ T , 


( 22 ) 
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can be expected to be noisy, bnt with (<h) = 0. 

As proper time increases the field rolls into a minimnm with = and dnring this 
“rolling phase” the pion mass can be imaginary, leading to nnstable growth of the Goldstone 
modes]^, |^. Since, as mentioned in the Introdnction, the energy density from the intrinsic 
chiral symmetry breaking is small^, |^], and the fireball snrface isolates the interior from 
the exterior of the light cone, it is reasonable to expect that the interior of the light cone 
ends np in a disoriented vacnnm. 

At late times snch a region of disoriented vacnnm with a given isospin orientation wonld 
relax back to ordinary vacnnm ($) = (a) = f.^, radiating pions with the same isospin 
orientation. 

It is also reasonable to expect that approximations based on the replacement of the fnll 
linear sigma-model by the simpler nonlinear sigma-model, with Lagrangian density given (in 
the chiral limit) by 




with cr^ -I- vf^ = 


(23) 


conld be reliably nsed at times late enongh so that the chiral field has already rolled into a 
minimnm with <1)1$ = (T^ -|- 7?^ = /^. 

In modeling these stages of evolution, the chiral limit can be safely taken as long as 
the proper time is small compared to m~^, while at proper times of order (1 to 2)m“^ the 
pion mass can no longer be neglected. At sufficiently large proper times, one should |]^ 
decompose the DCC field into physical-pion normal modes and let them propagate out to 
infinity as free states. 

Returning to the coherent state formalism reviewed in Section p.l| , we note that in the 
linear sigma-model we are dealing with four scalar fields. As a result, the appropriate “out” 
state is a coherent state characterized by four current densities 


|J«) = exp 


d~k 


q :=0 


exp 


i dk Ja{k)d\{k) 

. Q=0 


| 0 ) , 


(24) 


with the corresponding generating functional 


G[zq, Zi, Z 2 , Z:i] = exp 


dk Y.\Uk)Wzo.{k)-l) 


Q:=0 


(25) 


When not interested in observing a quanta (or working within the nonlinear sigma-model) 
one can simply set its fugacity variable zq to unity. This amounts to setting to zero the source 
associated with cr, which we choose to be Jq. This leaves the generating functional for the 
TT^, TT^ and TT^ fields 


G[zi^ Z2-, z^ 


exp 


dk j2\Uk)?{z,{k)-l) 

i=l 


(26) 


The physically observed fields are tt*^, 7r+ and tt , which are associated with the creation 
operators 
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Introducing the corresponding source currents^] 

Mk) = Uk) 


J±{k) = -^{Ji{k)TiJ2{k)) , 


one finds that 


'^Ji{k)a\{k) = Jo(fc)ao(fc) + J+(fc)a+(fc) + J_(/c)al(/c) , 


i=l 


and the generating functional for the tt’*', and vr fields is given by 
G[zo,2+,2-] = exp jdk (\Ja(k)\'‘(za(k) - 1 ) 

+\Mknzz(k) -1)+\j.(k)nz.{k) -1 


(29) 

(30) 


(31) 


(32) 

(33) 


For the DCC picture advocated in the following, the class of sources 


Ji{k) = \J{k)\ni , 


(34) 


where i = 1, 2, 3 and the n* are real constants such that 
From (|3^) one clearly has 


1, is of particular interest. 


Jo{k) = \J{k)\n 3 , 

J±{k) = l^^^(ni=Fm2), 
and the generating functional takes the form 


G[zo,z+,z_] = exp 


dk\J{k)\^i^{n3)\zo{k)-l) 


[{z^{k)-l) + {z.{k)-l)] 


(35) 


(36) 

(37) 


By realizing that / = (us)^, with / being the neutral fraction as in the Introduction, one 
can rewrite this generating functional as 


G[zo,z+,z-] 


exp 


ldk\j(k)f h{zo{k)-i) 

+ ^-^-^l(zz(k)-l) + (z.(k)-l)]\ 


(38) 

(39) 


This generating functional is appropriate for the description of events with initial con¬ 
ditions parametrized by a given /. However, for DCC production one must average over 

^This Jo should not be confused with the source term for cr, eliminated in the preceding discussion. 
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initial conditions with the appropriate weights P{f)- In such cases one can introduce 
the following type of generating functional 


G[zo, z+, Z-] = 


dfP(f)exp 

I ( 1 -/) 


[ dk |J(fc)P (^f{zo{k) - 1) 
[{z^{k) - 1) + (z.ik) - 1)] 


(40) 

(41) 


In particular, assuming that the initial rij are distributed with equal probability over the 
= l sphere, we get the characteristic form 


Pif) 


1 


(42) 


for the distribution of /. 

Finally, we observe that the above generating functionals with independent description 
of 7r~^ and vr” production can be turned into generating functionals for charged particles by 
hxing 


Zch = Z+ = z- . 

For example, for the generating functional (^TD and the P{f) of 


one finds 


Zch\ 


df 


exp 


/o 2 ^ 

+ (1 - f){zch{k) - 1)} 


dk \J{k)mf{zo{k)-l) 


(43) 


Thus far we have averaged over chiral orientations of the classical source, but not the 
overall shape of the source J. Eventually this problem must of course be faced. If the 
fluctuations about the mean, classical J are Gaussian, there exists a well developed formalism 
for dealing with them [OT. Indeed the DCC average we have performed is in fact also 
essentially carried out inKef. m by Andreev, Pliimer, and Weiner. 


In the case of DCC, the fluctuations will arguably go beyond the Gaussian approximation 
even when the DCC fluctuations are generated by a Gaussian distribution of initial condition 
parameters. This is a big problem, beyond the scope of this paper. We will comment on it 
again in the concluding section. 


4 Simple nonlinear sigma-model DCCs 

4.1 Setting up the problem in terms of an auxiliary field 

The picture of DCC evolution given in the sections 1 and 3 can be implemented (although, 
at least at early times, one does not obtain an accurate quantitative description of the 
physical system) within the nonlinear sigma-model in the framework of the set of classical 
solutions|l|, 1^, which have the form 
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( 44 ) 













where Vf is a constant but otherwise arbitrary unitary matrix]^, which orients $ along a 
direction in isospin corresponding to neutral fraction /, and 6 is such that 

□0 = 0. (45) 

(We are for the moment specializing our analysis to the chiral limit m^ = 0.) 

For Vf = 1 (which corresponds to / = 1) these solutions describe chiral helds with tti = 
7r2 = 0 and 


7r°= fT, sin 6 

”= 'Jfi- = /-cos» 


(46) 

(47) 


Our analysis will exploit the simplicity of Eq. (^5l), and will be based on the relation 
between source and particle flux discussed in the previous section. One aspect that perhaps 
deserves clarihcation is the relation between the pion flux and the solutions of the evolution 
equation for the held 9 in presence of a 0-source 

ne{x) = Je{x) . (48) 


For this, it is important to observe that when the source term Jq is well localized in space 
and time, as it must be for the consistency of our approach, the amplitude of the held 0 
must (because of energy conservation) become small everywhere at sufficiently late times. 
Therefore the relation between the late-time asymptotic helds and 0out, which according 


to (^^ is given by| 


= fn sm9out{x) , (49) 

can be well approximated by 

- U ^out{x) . (50) 

Based on the analysis reported in the previous section, it is then easy to see that the 
corresponding pion production is described by the generating functional 


G[z7r] 


exp 


dk f^\Je{k)\\z^{k)-l) , 


(51) 


and in particular, the inclusive spectrum of pions is given by 

dN^o fJ ~ , Vm o / N 

52 

d^k {2jTy 

For V/7^1, a simple generalization of this argument holds, and the (charged and neutral) 
pion production is ultimately described within the generating functional formalism discussed 
in the preceding section. 

§Note that the asymptotic behavior of the fields and sources plays a rather central role in our analysis. 
In particular, a class of solutions more general than given by $ = /tt with U ^ IT, could 

be considered, but would require everlasting sources. Such idealized sources have been considered in the 
analysis of Ref. ||T^ (where the source never turns off as a result of boost invariance) and Ref. Q| (where the 
source never turns off as a result of the infinite size of the “pancake”) and can be useful in deriving some 
intuition about DCCs, but do not reproduce the experimental conditions of DCC searches. 
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4.2 Pion flux for a class of sources 

We are finally ready to define in more precise terms a simple model which embodies the 
essence of the baked-alaska scenario. We start by introdncing a sonrce of the general formal 

Je{x) = A'K‘^f{t)Dret{.x) , (53) 

where 

Dret{x) = ^0(t)(5(r-t) (54) 

is the retarded Green fnnction for a massless scalar field, and f{t) is a fnnction of t 
such that 


fit) —0 t > cx) 

f{t)^l (55) 

i. e. the source is on at early times but eventually switches off. We are also using standard 
notations for r = \x\ and for the step function 0. 

The physically interesting quantity, as discussed in section |^, is Je{k), and the Fourier 
transform of the current density (^^ is given by 


Jg{k) = I dfx e^^^Je{x) 


, ifcOt 1 

at e — / dr r sm nr 


K Jo 


TT 


f{t)-Q{f)5{r -t) 


/-oo 


K Jo 

47^2 ^oo 


POO _ „ roo 

/ dt/(t)e**'* / dr 5{r — t) sin nr 
Jo Jo 


dt /(t)e*^'’*sinKt , 


K Jo 


(56) 


where we used the notation K=\k\. It then immediately follows that 


Mk) = 


dvr^ /-oo 


K Jo 


dt f{t)e 


iy/ 


sin nt 


(57) 


In particular, Je{k) can be approximated at large n as 


Je{k) 


27i^i /■'» , ^ 27i^iT 

2 k rsj - 


K Jo 


dt f{t)e^-^ 


K 


(58) 


where T = Jdt f{t). We see that the large k, behavior is not very sensitive to the exact form 
of the function f{t). However, this does assume that the sources are localized exactly on 
the light-cone. If we introduce a representation of the delta function with a finite width, 
the resulting field distribution is smoother and, as one might expect, the high frequency tail 
falls off much more rapidly. 

^The normalization factor has been chosen so that (as explicitly shown later on) for small enough t the 
DCC is pure tt = fn inside the light-cone. 
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4.3 Pion flux and fleld evolution for a specific source 

It is interesting to consider the idealized case f{t) = 0(T—t), where the sonrce term switches 
off snddenly at some time t = T, to be associated with the deconpling time. In this case the 
sonrce takes the form 

Je{x) = Q{T - t)-Q{t)6{r - t) , (59) 

r 


and Eq. (^) describing J off the mass shell rednces to 

Jeik) = - / dt sin nt 

K Jo 

27 r 2 r eiT{k°+K) _ I ^iTik°-K) 

K I + k) i{kP — k) J 


(60) 


The valne of J on the > 0 mass shell is easily obtained by snbstitnting + m? for 

and, in particnlar, in the massless case (and the large-K limit of the massive case) one 
finds 


Je{k) = -i 


A 2 T K. 


- 1 


2iK 


T 




K 

i ^{- 2 TK + i(l - 


(61) 


For the simple choice of sonrce considered in this snbsection it is also possible to derive 
explicitly the corresponding ^-field. We start by observing that 


9{x) = J d'^x' Dretix — x')Jg{x') 

J [ -1') 


= / d^x' 


,, [ 7r0(t') 


y^5{\x-x'\ -(t-t'))| 6{r'-t') 

1 


47r|x — X 

J dSQ{T-t')e{t')Q{t-t') (- 


x-y\^ 


^(|T— x'\ — (t — t')) 6{r' — t') 


= 'i*yQ(T - i')e(«')0(i - t') f S(p-(t- «')) S{r' - t') , 


(62) 


where in the last line we have introdnced the variable p = \x — x'\. We now go to spherical 
polar coordinates, choosing the z-axis to lie along x so that p = \/r‘^ + — 2rr'p, and find 

that 

6{x) =— dt'Q{t — t') / dr'6{r' — t') / dp—S{p — {t — t')) 

2 Jo Jo J-i p 

7 ]- pT fOO rr+r' 

= — / dt'Q{t — t') / dr'dir' — t') / dpSip — (t — t')) 

2r Jo Jo J\r-r'\ 

TV 

= — dt'Q{t-t') dr'6{r'-t')Q((t-t')-\r-r'\)Q{r + r'-it-t')) 

2r Jo Jo 

rp 

[ dt' Q(t — — t') — \r — t'\)Q{2t' — t + r) 

Jo 


TT 

2r 

71 


2r Jo 


[ dt' Q{t — t') {0(t — r)0(2t' — t + r) — 0(2t' — t — r)} 
Jo 


(63) 
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Here the last line (|63D follows after some straightforward manipulation of the step functions 
in the integrand. In order to render explicit the structure of the result (|63D it is convenient 
to examine separately the cases t<0, 0<t<T, and T<t] this also allows us to see how the 
various stages of the evolution of DCC described in the sections 1 and 2 are realized within 
this solution. From Eq. (|6^) it follows that 

• 9 = 0 when t<0, as expected since Dret{x) and Je{x) vanish for t<0. 

• For 0 < f < T 


e[t,r) = |0(t-r), (64) 

i.e. a uniform region of DCC is present inside the light-cone during the interval when 
the source is on. 


• Finally for times later than T 

TT / 27"' — / -4- 7^ \ 

e{t,r) = - ^j0(2T-t + r)0(t-r) 

(65, 

which is of the form [f{t — r) +g{t + ry\/r and therefore satishes the free wave equation, 
as expected since in this region the current density is off. The Eq. (^) also explicitly 
shows that at t = T, when the source is turned off, the ordinary vacuum (the a- 
vacuum) starts breaking into the interior of the light cone. At t = 2T the ordinary 
vacuum reaches the small-r region, and for times later than 2T the Eq. (|65l) can be 
interpreted as describing the outward propagation, in ordinary vacuum, of a localized 
DCC-wave. This sequence of events can be seen in the tt and a fields plotted in Fig. |I[ 


Having obtained the nonlinear sigma-model solution corresponding to the source (|5^), 
we can use it for a consistency check for the relation between source and held discusse^n 
section ^ i. e. we can check that using Eq. (^Of) one can indeed obtain the Fourier transform 
of the source (|5^) from the late-time 6*-held, which according to (|6^) is given (fot t > 2T) 

byQ 


^out 


TT /2T — t + r 
4 V r 


^ 0(2T-t + r)0(t-r) , 


( 66 ) 


We start by evaluating the three dimensional Fourier transform of dout 


dvr /■°° TT f2T — t + r 

— / dr r sm nr — 

4 


TT' 


K Jo 
2 


e{2T - t + r)e{t - r) 

roo 

/ dr (2T — f-|-r) sin Kr 0(2T — f-|-r)0(f — r) (67) 

Jo 


TT 


K Jo 

2 rt 


/ dr (2T — t + r) sin nr 

Jt-2T 


K Jt- 
2 


TT- 




^ sin^f sinK(2T —f) 

-2T cos Kt +-+-^^ 


K, 


K 


IINotice that the late-time 0-field satisfies \0out{t,r)\ < •nT/{2t), and therefore the approximation TTout = 
A sin dout ~ f-nQout holds asymptotically, as argued at the beginning of this section based on general energy- 
conservation arguments. 
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The three dimensional Fourier transform of 6 is then given by 




dt 


k) 



2T sin K,t H—(cos nt 

K 


cos k(2T — t)) 


( 68 ) 


Finally, following Eq. (|^) we get (also taking into account that this calculation is in the 
massless limit, and therefore ujk = f^) 


Jeik) = —i e 


int 




(3) 


= —I e 


IKt 


TT 




2 r 






71 


—2Tk + i 


qi- 


K 

^i2TK,\ 


K 


(69) 


in complete agreement with Eq. _ 

In Fig. 1 we show a few snapshots of the evolution of the helds for the solution (pS]) and 
the corresponding source obtained using Eq. (Ell), which agrees very well with the expression 

( 0 )- 


4.4 Correlation between generic and DCC pion production 

The question of correlations between generic and DCC pion production is very important 
for experimental DCC searches. It is quite reasonable that such correlations should exist. 
In our baked-alaska scenario, if the source strength on the light cone is large, i.e. lasts for a 
long time T, then the amount of DCC which is produced will, as we have seen, also be large. 
To be more quantitative is not easy. Here we attempt to make the connection by a simple 
argument based on energy densities. In the absence of the source, then there would be a 
surface-tension, i.e. an energy per unit area, associated with the boundary region between 
the sphere containing the DCC and the normal vacuum on the outside. It is contributed by 
the kinetic-energy term of the DCC hamiltonian, since gradients of the pion field exist in the 
interface region. It is reasonable to assume that when the energy per unit area of the generic 
partons or hadrons in the source region exceeds this surface energy, the inner DCC will be 
decoupled from the outer vacuum, and conversely when this energy is less than the DCC 
surface-tension, the source term will be inoperative.^^ This will allow a connection between 
the amount of generic production and the amount oiDCC produced. 

In the spherically symmetric case, by assumption the DCC surface energy is contained 
in a shell of thickness A, with A a typical hadronic scale, say 0.2/m < A < 0.6/m. We can 
therefore write 


Edcc 

A 




dvrT^ JT- 


(70) 


**It may also be argued |0| that the source of the pion field should be a scalar density built from the 
constituent quarks composing the generic material on the light-cone. This quantity is not the energy- 
momentum tensor, so that this is not obviously the same criterion as we are using. On the other hand, when 
the fields are rapidly varying (as they are in the source region), it is not clear what the correct choice of 
source of pion degrees of freedom is, and a simple argument based on energetics seems not unreasonable to 
try. 
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Figure 1: Evolution of the vr (continuous line) and a (broken line) fields according to the 
classical nonlinear sigma-model, starting from pure DCC tt = f-K inside the light-cone at 
t = 1/m. The helds are shown in the snapshots on the left, whereas the snapshots on 
the right show the modulus squared of the pion source function (the J of Eq. (p0|)). The 
horizontal scales for the left-side and right-side plots are in fermi and (fermi)“^ respectively. 
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where the last approximation on the right-hand side follows from approximating the prohle 
at decoupling of the 6*-£eld near r = T with a linear interpolation between 9{r = T — A,t = 
T) = 1 and 9{r = T,t = T) = 0. 

For the energy per unit area of the hot shell made of collision debris, one easily hnds 


Eghell dEsfiell 

[<PT>^dni 


< Pt > dN 

A dA 

T^dn 

generic 

T2 dQ_ 


(71) 


generic 


If indeed at t = T the DCC surface ener gy d ensity equals the energy per unit area of the 
hot shell of collision debris, one finds from (ffoD and (Q) that 


<Pt> 


dN 

dVt 


= T 


generic 


2 J-K 

2A 


(72) 


The aforementioned correlation between generic and DCC pion production is then seen 
upon combining ([7^) with (|^ and Specifically, one finds (in the chiral limit) 


P 


dN 

dQdp 


J DCC 


2ti 

'a 


<Pt> 


dN' 

dVt 


(73) 


generic 


We therefore see that it is possible for DCC production and generic production to be 
comparable in terms of the number density of produced particles. 


5 Beyond the classical nonlinear sigma-model 


Up to this point we have based our description of DCCs on linear, semiclassical, coherent- 
state solutions of a simple nonlinear sigma-model. However, there has been a considerable 
body of work on DCC production which goes substantially further. Not only are the nonlinear 
equations of the linear sigma-model (or even more complicated models [Q) considered, but 
also the effects of quantum fluctuations are taken into account within the mean-held [O, ^ 
(or large-iV O, |H|) approximation. This level of calculation has become the de facto ^^ate- 
of-the-art”. However, in most cases the space-time geometry is greatly simplihed, or else the 
approach is aggressively numerical. 

The closest calculation at this level to what we have presented here has been performed 
by Lampert, Dawson and Cooper (LDC) |^. They consider a boost-invariant spherical 
expansion, such that the helds depend only upon the proper time which has elapsed since 
the expansion began. This is not very realistic, because the inclusive particle distribution 
which emerges must be the same in all reference frames, and therefore requires an inhnite 
mean energy per particle, and an inhnite formation-time for the hnal-state distribution. 

While the LDC solutions are, as they stand, not very realistic, it is not too hard to adapt 
them to the baked-alaska scenario which we have described. We shall sketch in subsequent 
subsections how this works. The main point is that if we assume that the dynamics is as 
described in LDC for times t less than T, after which time the sources on and near the light 
cone are turned oh, then by causality the LDC solution will still be exact within the double¬ 
cone region, i.e. between the forward light-cone (vertex at t = 0) and an inverted light cone 
with vertex att = 2T. If T is large enough, one could hope that the helds inside the double¬ 
cone region and far from the light-cone be asymptotic, so that they could be matched onto 
the nonlinear sigma-model helds we use, and an estimate of the low-momentum portion of 
the particle spectrum could then be obtained using Eq. (^H). The quality of this method can 
be tested by varying T and determining which portion of the inclusive spectrum is insensitive 
to T. 
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For a more rigorous analysis of the spectrum, the source should be turned off at the 
decoupling time T, but the helds should be evolved according to the linear sigma-model 
up to times late enough for the evolution to be effectively free. At sufficiently late times 
one can then reliably extract the information on the inclusive spectrum using the procednre 
described in Sec. 4. 

It is also instructive to consider the simplification achieved by neglecting the mean-field 
qnantnm corrections, in which case the LDC calculation is rednced to the solntion of conpled 
ordinary differential eqnations describing the evolntion of the classical fields in proper time. 
We have made snch calcnlations for initial conditions chosen by LDC, and find remarkably 
close agreement of the time dependence of the pi and sigma fields with what is obtained 
from the fnll mean-field quantum calculation. This is enconragement that, when one goes on 
to consider the linear sigma-model in more difficult, less symmetric geometries, a classical 
calculation may well snffice to provide at least a semiqnantitative pictnre of the dynamics 
which the mean-field calculation would provide. After all, only a semiqnantitative description 
need be obtained from the linear sigma-model, because it is just a rough approximation to 
the complete low-energy effective chiral Lagrangian of real QCD. 

In the following snbsections we sketch more details of this line of argnment. In snbsection 
5.1 we discuss the connection between the pion flux as calculated from the nonlinear sigma- 
model with what one would obtain from a classical solution of the linear sigma-model. 
In subsection 5.2 we describe in more detail the LDC analysis, especially in the classical 
approximation, and we nnmerically compare the mean-field and classical LDC-type solntions. 
In subsection 5.3 we match the classical versions of the LDC solntions at time T onto the free 
asymptotic fields of the nonlinear sigma-model, thereby defining an effective source function 
J, from which the pion distribntions are calculated. Finally, in snbsection 5.4 we match the 
LDC solntions at time T onto fields of the linear sigma-model, evolve according to the linear 
sigma-model np to some time T' say (late enongh for the evolntion to be effectively that of 
a free field) when the calculated effective source function will represent the actnal pion flux. 

5.1 Deriving the pion flux in more general frameworks 

We start by applying onr formalism to the derivation of the pion flnx for the fnll linear 
sigma-model or related models, rather than for the nonlinear sigma-model considered in 
the previous sections. Provided one is dealing with a well-defined scattering problem, with 
the sonrces localized in space and time, the evolntion of the n and a degrees of freedom 
will eventnally reach an asymptotic regime governed by free field behavior. An estimate 
of the associated inclusive pion spectrnm can indeed be obtained by applying the formulas 
discussed in the previous sections. 

While technically this procednre is rather straightforward, it is important to realize that 
the associated “sources” are somewhat different from the ones we have been discussing. In 
this more general case one is actnally dealing with “effective sources”, useful as computational 
tools in the analysis of pion prodnction, bnt not to be interpreted as physical external sources 
in the problem. From the point of view of the original model, say the linear sigma-model, 
these effective sonrces are given by the sum of a physical external source and a term from 
the self-interactions of the fields. 

These ideas are of rather general applicability; for example, in the investigation of an 
interacting system described by the Lagrangian density 

1 HI ^ 

C = - —$2 - D(<F) (74) 

one is natnrally led to the stndy of the evolntion eqnation 

(□+ m2)4> =-W($) + J , (75) 
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where J is a physical “external” (^-independent) source and J^ff = —V'{^) + J is an effective 
source. 

The simulation reported in Fig. corresponds to the linear sigma-model classical evolu¬ 
tion from an initial configuration given by pure DCC, n = inside the light-cone and true 
vacuum outside (he. a snapshot of the solution (^) at some chosen time) and vanishing 
initial field velocities (except on the light cone). We simulate the classical field equations, as 
obtained from the Lagrangian density (^I]), for a spherically symmetric field configuration. 
Rather than simulating the <h field directly, our program evolves r$(f, r) which simplifies 
the form of the d’Alambertian. The boundary conditions at the origin are set up to ensure 
that <F(f,r) is an even function of r at all times. The relevant classical field configuration is 
then evolved in time, using a simple staggered leapfrog algorithm (see for example |Q). The 

Fourier transforms involved in finding the source current J{k) as defined in expression (pOD 
are done using a straightforward (extended trapezoid) method]^ Fig. 0 shows a selection of 
output snapshots chosen to illustrate the main features of the evolution. The pictures on 
the left describe the evolution of the pion and sigma fields, while the pictures on the right 
describe the corresponding “evolution” of the Fourier-space effective source function. The 
emergence of a stationary Fourier-space effective source (which encodes the information on 
particle production) reflects the fact that at late times the evolution of the tt and a degrees 
of freedom reaches an asymptotic regime ruled by the nonlinear sigma-model. However, it 
should be noted that the low momentum part of the spectrum is only complete after a time 
scale between 4 and 8 fermi, and is significantly larger (a factor of 2) than what was obtained 
for the nonlinear sigma-model for identical initial conditions. 


5.2 Classical version of LDC approach and reliability of coherent- 
state descriptions 


Our description of DCCs uses the classical equations of motion to obtain an “out” field 
from a given “in” field. This “out” field is then mapped into a corresponding coherent 
state from which particle production (a quantum effect) can be derived. Some elegant 
recent studies 0.00 have been based on more general formalisms for the description 
of quantum effects and~have taken into account (some of) the non-perturbative quantum 
effects contributing to the structure of the full propagator. One is then confronted with the 
solution of a genuinely non-classical evolution problem, in which “gap equations” describing 
the full propagators are combined with (modified) evolution equations for the fields. Seen as 
solutions of a variational problem, these equations result from finding an extremum of the 
(quantum) effective action for composites |p5[|, just like the classical evolution equations are 
obtained from finding an extremum of the classical action for (local/non-composite) fields. 

In Ref. 1^ , LDC investigated the chiral phase transition by modeling the relevant hadron 
dynamics with a linear sigma-model, and adopting evolution equations that take into ac¬ 
count part of the non-perturbative quantum effects contributing to the structure of the full 
propagator via the familiar large-formalism. They concentrated on boost-invariant spher¬ 
ical expansions, such that the mean-field expectation values depend only upon the proper 
time T = — r^. We refer the interested reader to Ref.for the complete description 

of the LDC approach. For the purposes of the analysis presented in the remainder of this 
section, it is sufficient for us to consider explicitly the evolution equations corresponding 
to the classical version of the LDC equations, i.e. obtained from the LDC equations by 


^^Note that by setting X = H = 0 in (|2^) we obtain the free wave equation. This enables us to use an 
almost identical program to evolve the 0-field in the nonlinear sigma model to obtain the graphs in Fig. |^. 
The major change in the code between the two cases is in fact the use of the identity tt = /tt sin 9 in the 
J-extracting routine. 
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dropping all contributions coming from the dressing of the propagator: 


1 


dr 

1 dr) 


'1 d 


dr 

1 dr) 


o- = f^ml 


TT = 0 


(76) 

(77) 


As a preliminary test of the reliability of our description of DCCs based on the classical 
evolution equations, we have compared the results for the mean-held evolution of the tt and a 
helds obtained in Ref.|]l8|, to the corresponding results from the classical equations ([76D-([r?|). 
In Fig. 3 the results of this comparative analysis are reported for the initial conditions sin^d 
out in Ref. ilS 


as the most “DCC-favoring” within the special family of initial conditions 


considered there; specihcally, we integrate the helds from the initial conditions 


a = 0 fm ^ & = —1 fm ^ 

TT = 0.3 fm~^ TT = 0 fm~‘^ (78) 


at r = 1 fm. Since the focus of this exercise is only on the held evolution, rather than particle 
production, for simplicity we kept the exact (spherically-symmetric and boost-invariant) 
source structure adopted in Ref.[^. Fig. 3 suggest that even at a quantitative level the 
description might be satisfactorily accurate. This is especially so because one is in any 
case using rough models of the relevant hadron dynamics (i.e. it appears to be likely that 
the inaccuracies introduced by using classical evolution equations might be less important 
than the ones resulting from modeling the relevant hadron dynamics with, say, the linear 
sigma-model). 


5.3 Low-momentum portion of the inclusive spectrum in the LDC 
approach 

As explained at the beginning of this section, the spherical expansion investigated in Ref.||18||, 
which is fueled by everlasting sources, and involves helds depending upon only the proper 
time, is not very realistic. Still, as mentioned above one could attempt to extract the 
low-momentum portion of the inclusive pion spectrum associated with a baked-alaska-type 
modihcation of the LDC approach, by mapping the LDC solutions at time T onto free 
asymptotic helds of the nonlinear sigma-model, thereby dehning an ehective source function 
J, from which the pion distributions are calculated. Ideally, one might hnd that for large 
enough T the helds inside the double cone region be almost everywhere asymptotic, and 
that only the high-momentum tail of the particle spectrum could not be captured by such 
an approach. In practice, however, we hnd that not even the portion of the inclusive pion 
spectrum with very low momentum is well determined at times as late as 50 fm. 

In Fig. 4 we report the results of one such analysis, in which we integrate the helds from 
the initial conditions (^) at r = 1 /m, to some later time T and use the boost symmetry 

to reconstruct the helds everywhere on the t = T surface inside the light cone. J{k) is 
then extracted from this held conhguration using Eq. (pOD as before. The a and vr held 
conhgurations and the ehective source function described for various choices of the above- 
mentioned time T in Fig. 4 clearly rehect the shortcomings of the approach discussed in this 
subsection. 
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Figure 3: Evolution of the pion and sigma fields as functions of proper time in the LDC 
setup. The continuous line corresponds to the purely classical analysis, whereas the dotted 
line corresponds to the quantum analysis reported in Ref. [18]. At proper time r= 1/m the 
helds and their derivatives are fixed to be tt = 0.3/m“/ a = 0fm~^, n = 0fm~‘^, & = 
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Figure 4: Evolution of the tt (continuous line) and a (broken line) fields, and the pion source 
function according to the classical linear sigma-model, in the presence of a boost-invariant 
and spherically invariant source considered in Ref. [18]. Again, the initial conditions are 
chosen according to the “DCC-favoring” scenario considered in Ref. [18]. 
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5.4 LDC approach with truncated sources 


The deficiencies of the method discussed in the previous section can be easily remedied. 
Evidently, after the external source is turned off at the decoupling time T, the fields should 
be evolved according to the linear sigma-model up to times late enough for the evolution to 
be effectively ruled by free field behavior. At such late times one can reliably extract the 
information on the inclusive spectrum using the procedure described in Sec. 4. 

In Fig. 5 we report the results of such a simulation, again as snapshots describing the 
evolution of the fields and the effective source function corresponding to the initial configu¬ 
ration singled out as 


‘DCC-favoring” in Ref. 
simulation a large decoupling time (5/m). 


For illustrative purposes we chose in this 
IrT^this case the effective source does reach a 
stationary regime; however, by comparison with Fig. 2 we see that this asymptotic behavior 
only emerges at rather late times (~ 20 — 30 fm). The comparison of Fig. 4 and 5 shows 
that the approach discussed in the preceding subsection, in which the LDC sources were 
never turned off, can largely overestimate {e.g. by more than a factor 10 if T < 5/m, as 
assumed in Fig. 5) even the low-momentum portion of the spectrum. 


6 Conclusions 


In this paper we have investigated the various stages of the evolution of disoriented chiral 
condensates via the “baked-alaska” mechanism. Most of our analysis has been elaborated 
using classical equations of motion based on either the linear or nonlinear sigma model. The 
associated framework of coherent states was then used to make the connection with the dis¬ 
tribution functions for the particle production. Important in this step is the identification of 
the source-function of the produced particles, namely the right-hand side of the usual wave 
equation (cf. Eq. (|75D ). The square of the on-shell fourier transform of this source func¬ 
tion, as determined from the solutions of the equations, then provides directly the inclusive 
distribution of particles. 

In general the source term consists of two parts. One is concentrated near the light cone, 
and is a genuine external source, not a function of the chiral fields, to be associated with 
the generic collision debris of partons, constituent quarks, etc. The other part consists of 
the nonlinear terms, built from the chiral fields themselves, which appear in the classical 
wave equation. We found evidence, within the nonlinear sigma-model approach, that the 
number of produced DCC pions is likely correlated with the number of generic hadrons 
produced, with this correlation local in (lego) phase-space. The number of DCC pions could 
be comparable with the number of generic hadrons according to this crude estimate, but the 
uncertainties are very large. 

We also compared our very simple classical approach with a mean-field calculation which 
includes one class of quantum corrections, and at least in the case we studied, the quantum 
effects appear not to be of great importance. This is encouragement that, when attempts to 
go beyond the spherical symmetry assumed in this work, the simpler classical approach may 
suffice to reveal most of the important physics. 

In all of the work in this paper (and in most of the literature), we have assumed spherical 
symmetry of the solutions. Regrettably, this geometry is too simple for many realistic 
applications. The intrinsic sources are reasonably uniform in lego variables, not spherical 
coordinates, and this geometry needs more detailed study. In addition, fluctuations about 
the mean behavior are very important. A piece of DCC with relatively large transverse 
velocity will look in the laboratory like a coreless minijet, with contents containing small 
relative momenta. So the source distributions most relevant to DCC searches in high energy 
hadron collisions should not only be described in lego variables, but also contain minijet 
clusterings. 

However, in defense of what we do, each piece of DCC in momentum space is a cluster 
of pions of near identical momenta — a “snowball” — which has a local rest frame [0, ITO . 
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Figure 5: Evolution of the tt (continuous line) and a (broken line) fields, and the pion source 
function according to the linear sigma-model, in the presence of a trnncated version (as 
described in the text) of the a boost-invariant and spherically invariant source considered in 
Ref. [18]. The sonrce is switched off at f= 5/m. 
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In a snowball rest-frame, the calculations we make should be a reasonable description of the 
dynamics of that particular snowball. But one needs to know how the chiral orientations 
and probability of occurrence of snowballs which are neighboring in momentum-space are 
correlated. Very little work on this exists. 

In addition, one should average over sources more broadly. This includes not only the 
properties of the intrinsic sources discussed above, but also the initial conditions imposed 
on the chiral helds at early proper time, i.e. at the onset of the chiral symmetry breaking. 
A good starting point will be to do this for the classical version of the interesting model of 
Lamport, Dawson, and Cooper |^, truncated at large times as described in this paper. 

The hnal product of all this should be a generating functional for DCC particle pro¬ 
duction, which is an average over sources of a Gaussian generating functional characteristic 
of a coherent-state and classical solution produced by a specihc source (see, for example, 
reference for the formalism). However, even when this is attained, it still leaves open 


the more difficult problem of synthesizing such a generating functional with one for generic 
production, since there is not yet a consensus on what represents a good choice for the latter. 

The emphasis we make in this paper on DCC sources makes the formalism look more and 
more similar to what is used to describe Bose-Einstein correlations [^. There is certainly 
a close relationship What we believe special about the baked-alaska scenario, even ac¬ 


cepting only the broad outlines, is that there is assumed to be nontrivial dynamics occurring 
deep within the future light cone. Irrespective of details of our modeling, the presence of 
such dynamics would be a new element in the description of hadron-hadron collisions at high 
energies. 
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